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Abstract 

The interacting and holographic dark energy models involve two important quantities. One is 
the characteristic size of the holographic bound and the other is the coupling term of the interaction 
between dark energy and dark matter. Rather than fixing either of them, we present a detailed 
study of theoretical relationships among these quantities and cosmological parameters as well as 
observational constraints in a general formalism. In particular, we argue that the ratio of dark 
matter to dark energy density depends on the choice of these two quantities, thus providing a 
mechanism to change the evolution history of the ratio from that in standard cosmology such 
that the coincidence problem may be solved. We investigate this problem in detail and construct 
explicit models to demonstrate that it may be alleviated provided that the interacting term and 
the characteristic size of holographic bound are appropriately specified. Furthermore, these models 
are well fitted with the current observation at least in the low red-shift region. 
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I. INTRODUCTION 



Recent astronomical observations indicate that our universe is currently undergoing an 



epoch of accelerated expansion 



Following the standard Friedmann- 



Robert son- Walker (FRW) cosmology such an expansion implies the existence of a dark 
energy (DE) component to the mass-energy density of the Universe'^. At present it is fair 
to say that disclosing the nature of DE is one of the central problems in the research of 
both cosmology and theoretical physics (for recent reviews, see Q^lsiyi)- In this direction 
we are faced with many fundamental and difficult issues, among of which the following 
two open questions are of particular importance. The first one is on the nature and dy- 
namical properties of the dark energy. Currently it is not clear yet whether DE can be 
described by cosmological constant which is independent of time, or by dynamical fields 
such as quintessence, K-essence, tachyon fields or phantom fields. The second is the coin- 
cidence problem, dubbed as " why are the densities of matter and dark energy of precisely 
the same order today?" lfl| . 

To shed light on these two open questions, some interesting DE models were proposed 
recently. Those models can be divided into two categories, i.e. the holographic dark energy 



(HDE) models and t 



hypothesis 



16 



le interacting DE models. The former stems from the holographic 



17| and can provide an intriguing way to interpret the dynamics of DE, 



while it is suggested that the latter can help to understand the coincidence pro 
considering the possible interaction between dark energy and cold dark matter 



21 



)roDiem by 



3- 

Let us first start with a close look on the holographic dark energy model motivated from 
the holographic hypothesis, which has gradually been believed to be a fundamental principle 
in the quantum theory of gravity. According to this principle, the number of degrees of 
freedom for a system within a finite region should be finite and is bounded roughly by the 
area of its boundary. While in a cosmological setting the challenge is to put a reasonable and 
well-defined upper bound on the entropy of the universe. Motivated by Bekenstein entropy 
bound, it seems plausible to require that for an effective quantum field theory in a box of 



^ For recent discussion on the possibility of constructing accelerating universe without dark energy, see for 
instance EH El HQ. 
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size L with UV cutoff A, the total entropy should satisfy the relation 

I^^^A' < Sbh = ttL^MI (1) 

where Sbh is the entropy of a black hole with the same size L. but further consideration 
indicates that to saturate this inequality some states with Schwarzschild radius much larger 
than the box size have to be counted in. As a result a stronger entropy bound has been 
proposed in jl^ , requiring that the total energy of system with size L should not exceed the 
mass of a black hole with the same radius, namely 

L^A^ = LVa < LMl (2) 

While saturating this inequality by choosing the largest L it gives rise to a holographic 
energy density 

PA = 3c'M'^L-', (3) 

where c is a dimensionless constant. Then the key issue is what possible physical scale one 
can choose as the cutoff Lconstrained by the fact of the current acceleration of the universe. 
Originally the natural choice is to identify the Hubble horizon as L, however, as pointed out 
in llGll, this will lead to a wrong equation of state for dark matter which conflicts with the 

n 

ordinary one in standard cosmology. As a result, in Li proposed to take the future event 
horizon as the largest size L, which gives rise to desired results and then stimulate a lot of 
interests and discussions in this subject However, there are some unsatisfactory points 
in this conjecture. First, it still remains puzzling how the current evolution of dark energy 
density can be determined by the future event horizon. Second, the coincidence problem 
can hardly be solved in this context. 

As pointed out in [2^, the reason that one is forced to take the future event horizon is 
based on the assumption that the energy densities of dark energy and dark matter evolve 
independently. However, if there exists interaction between DE and dark matter (DM), then 
the cutoff L is not necessarily identified as the future event horizon. As a matter of fact, the 
interaction between DE and DM is proposed to solve the coincidence problem and has been 



discussed in many recent works 2J| . This can be accomplished by introducing the coupling 
terms in the equations of state for matter and dark energy densities, which can bring the 
ratio of these two ingredients into a constant at late times. From the theoretical point of 
view this sort of coupling is completely possible due to the unknown nature of DM and DE. 



In addition, this proposal is compatible with the current observations such a s t 
WMAP dataQ, 

and even favored in some circumstances as suKKested in 



le SNIa and 
But until 



now, only certain special interacting terms have been considered in existing literatures. 

Now based on the above discussion, it is natural to ask if we could combine these two 
theoretical proposals together so as to improve our understanding of dark energy. This 
is the main purpose of our paper. Although there are many existing works in both two 
directions, most discussions in those works only considered specific characteristic sizes or 
interacting terms. For instance, the characteristic size is usually assumed to be the future 
event horizon after the work and the interaction term is assumed to take a form as Sb^Hp 
where 6 is a coupling constant. However, there is no strong theoretical motivations for these 
choices. In this paper, rather than fixing either the interacting term or the holographic 
characteristic size, we intend to investigate the nature of interacting and holographic DE 
and the coincidence problem in a more general formalism. 

Our paper is organized as follows. We first present brief reviews on interacting dark 
energy and holographic dark energy in Section 2 and Section 3, respectively. In particular, 
given the conditions that our universe is currently accelerating and the ratio of dark matter 
density to dark energy density deceases, we derive general constraints on the relations among 
the interacting term, holographic size and the equation of state parameter of dark energy u. 
Then we turn to the coincidence problem in Section 4, under the simplest requirement that 
the ratio of dark matter to dark energy density to be constant. In Section 5 we consider the 
case that the ratio can vary with time slowly and demonstrate how the coincidence problem 
can be alleviated through some specific examples by providing appropriate interacting terms 
and holographic sizes. 



II. INTERACTING DARK ENERGY 



We start with the standard Friedmann equations in which DE and DM are assumed to 
be independent and there is no interaction between them. Provided our current universe 
is dominated by dark energy with the state equations p = up^ and cold dark matter with 
p = 0, these equations read as 

PA + PM = ^MlH\ (4) 
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H = --M~^[Pa{1 + uj) + pm], 



(5) 



Pa + Pm = -3H[p a{1 +u;) + Pm], (6) 
where H = a/a is the Hubble factor. It is well known that equations (jH)-® are not 
independent and any one of them can be derived from the other two. Introduce f^A = 
Pa/ (3Mpif^) and Qm = PM/(3Mpi/^), the first Friedmann can also be written as ^a + ^m = 
1. 

Now we proceed to interacting dark energy models in which dark matter and dark energy 
are postulated to be coupled such that dark energy can decay into cold dark matter. As 
a result, the last equation can be written as the combination of following two evolving 
equations, 

Pa = -3Hpa{1 + uj)-Q, (7) 



Pm = -3HpM + Q, (8) 

where Q denotes the interacting term. To be a realistic model, the interacting DE model 
should satisfy the observational constraints. First, we consider the constraint on Q by the 
observation that our universe is currently accelerating. Since the ratio of dark matter to dark 
energy plays a special role and its dynamics is a major subject of this paper, for convenience 
we denote the ratio Pm/pa by r which is related to Qa by 1 + r = 1/Qa- Then from (j^ and 
(0) we have 

Noticed that equation © always holds no matter whether the interaction is taken into 
account or not. The solution to this equation can be formally written as 

i/ = i/oe"^^"'('+^)"^ (10) 

where x = Ina. As a result, the requirement a > leads to 

l+r + 3cj<0. (11) 

On the other hand, from (|7j) and (jH)) we find that the interacting term has the following 
general form. 



where H can be absorbed by redefining r / H = dr / dx = r' ^ such that 

Q = Y^(r'-3c^r). (13) 

Then from one finds that 

Q>r + -^. (14) 

1 + r 

Furthermore, it is expected that the ratio of dark matter density to dark energy density 
decreases with the evolution of the universe, namely, r' < 0. This requirement together 
with the previous one (jllj) imphes that the interacting term should satisfy the following 
constraint 

^' A -Scur . ^, 

r + - <Q<- . 15 

1+r 1+r ^ ^ 



III. HOLOGRAPHIC DARK ENERGY 



Now we turn to the holographic dark energy (HDE) models. As introduced in Section 1, 
in this context the dark energy density is assumed to be saturated in the region with size 
L, 

PA = Sc'M'^L-'. (16) 

Comparing this bound with the first Friedmann equation we easily obtain a relation between 
the characteristic size L and the Hubble factor H as 

LH = VTTrc. (17) 

In general the characteristic size L need be chosen in such a way that the dark energy 
can be responsible for the acceleration of the universe. Furthermore, such choice should 
not conflict with the state equations of dark energy and dark matter. As a result, when 
there is no interaction between DE and DM, i.e. Q = 0, L is conventionally taken as the 
future event horizon 0, |3] so as to fit the observational data. However, as explained 
in Section 5, the coincidence problem is hardly solved by pure HDE alone. Therefore, in 
the following discussion on HDE models, we will consider the case with non- vanishing Q. 
Moreover, instead of fixing L or the interacting term Q at the beginning, we will take a more 
phenomenological view and consider them to be free dynamical quantities but constrained 
by observations. 



6 



Now, as in the previous section, from the requirements of accelerating universe and de- 
creasing ratio of DM density to DE density one can derive the constraints on those quantities. 
Taking the derivative with respect to time on both sides of equation (fTTj) and using (fT^ 
lead to a relation between L and Q as 

Q = r'-2r{^-^-). (18) 

Now due to the constraint (jl5|) we find the size L should satisfy 

?) ,r' uj . L' r' , , 

2(3;: + TT7 + 1' ■= T ■= 2(TT7) + 1- ^^"^ 

Thus we find in a general formalism of interacting holographic dark energy, the interacting 
term and the characteristic size of holographic bound are constrained by the inequalities 
(fT3j) and (fT^. respectively. 

Alternatively from and ()18|) we may have a relation between L and u as 

2^ = 3(1 + + (20) 

Thus r' < leads to 

— <-l + . 21) 

Furthermore replacing the parameter uj appearing in (|T5|l by L we may find the following 
inequality 

"^^^L "2^^^^^^T + 2"^^- ^^^^ 
In summary, we find for interacting holographic dark energy models, L, Q, uj and r are 
not independent quantities but related by equations (fTSj) and (pUj). or (fT^ and (pUj) since 
among these three equations only two of them are independent. Thus given any two of 
them then the dyanamics of the other two can be determined. For instance, if we specify 
the interaction term Q and the characteristic size L, then the dynamics of u and r may be 
determined and vice versa. However if only one of them is specified, then the dynamics of 
the other three can not be uniquely fixed. For example, in an interacting HDE model with 
the future event horizon as the characteristic scale L, i.e. 
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which gives rise to a relation L = LH — 1. Then from ()20j) it is easy to derive the following 

relation 

1 \/n^ r' 

-(l+r + 3a;) + ^^ + -=0. (24) 
In this case r' < requires that 



l+^ + 3^>_^^/I±Z. (25) 

If we further specify the interaction term, e.g. Q = 36^(1 + r) as in references [2^. then the 
dynamics of uj and r can be determined uniquely, as discussed in j^]. However, since our 
goal is to investigate the coincidence problem, we intend to put constraints on the evolution 
of r and then explore what expressions the other quantities including L,Q and u may take. 
This is what we are going to do in next two sections. 

IV. COINCIDENCE 

Before proceed, we first demonstrate how the coincidence problem arises in the standard 
cosmology. Setting Q = in ((Zj) and (jH)) will lead to r' = 3ujr (see also ()13j) ). Now from 
f|TT|l one finds that 

dlnr 

— — = 3cj < -r - 1 < -1, 
ax 

which means during acceleration r decreases faster than a~^. Furthermore, from Friedmann 
equations, one finds that 

r = ro (a/ao)~^ , (26) 

and 

a/ao = C(e^*/*«-e-^*/*«)'/', (27) 

where C and A are 0(1) constants which can be related to the current values of Q/<^. Then 
it is easy to see when t Iq, r oc and r decreases quadratically as expected for a matter 
dominated universe and when t ^ Iq, r deceases exponentially as expected in a dark energy 
dominated universe. Then it is only when t is around to that r ~ 0(1). 

It is expected that adding interaction may change the dynamics of r greatly. In this 
section we consider the simplest possibihty with r = 0, which implies p\ oc pm oc H"^. It is 
worthwhile to stress that this situation only occurs at late times. Suppose the ratio r is a 
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constant, i.e. r = tq. We immediately obtain the following equations 

~ _ -3^ro 
1 + ro 

L = + (29) 

In addition, the Hubble factor is inversely proportional to characteristic size L as H = 
y/1 + rQcL~^. Therefore specifying any one of quantities Q,L,uj, the dynamics of the other 
two can be uniquely determined from above equations. We classify some possibilities in the 
following subsections. 



(28) 



A. r = with only interaction term specified 

One possible choice for the interaction term is setting Q = 36^(1 + ro) as in previous 
references, where is a constant. Then from (OHl) we find that u is fixed as 



, . .,/l±l2}l^ (30) 



ro 

Consequently the solutions to H and L can be obtained from as 



4(l_fe2_6^) 

1^ rn ' 



H = Hoa"^''-" --o\ (31) 
In addition, from ()15p one finds that 

36' > = ^M- (32) 

(1 + ro) 

Therefore, it is obvious that in non-interacting models (i.e. 6 = 0) r = and a > cannot 
be achieved simultaneously. This can be considered as an important hint for the need of 
interacting dark energy. 



B. r = with only holographic characteristic size specified 

As shown in the previous section, specifying the holographic characteristic size will de- 
termine Q and u since r has already been fixed to be tq. Here we consider the HDE model 
with L being the future event horizon. From (f^H|) and ()29|1 we find correspondingly that the 
parameter u and Q have to be constants as well. 

Q = Qo = roil+ J— )■ (33) 
9 



c. = -o = -^(l + ro + ^^^). (34) 
Using the current data pmo — 0.25, pao — 0.72 and setting c = 1, we find 

ro - 0.35, cuo - -1.22, (35) 

which is a phantom-preferred model. From we find that 

2(l+'-0) , , 

a ~ t3{l + ro+a,o)^ (36) 

while 

Pa ~ Pa/ ~ t'"^- (37) 

C. r = with both interaction term and holographic characteristic size specified 

If both L and Q are specified as previous subsections, then from pn|l and (jH^ it is easy 
to see r = can be reached only when the constant 6^ takes the value as 

'2 ""^o^o - 0.24. (38) 



V. SOFT COINCIDENCE 

The above discussion shows that the interaction between DE and DM can lead to a 
constant r. Although it is not clear how to obtain a r of 0(1) size at early times, this simple 
strategy can be used to account for particular situations (e.g. late time evolution of the 
universe). Nevertheless, there is no strong motivation for setting r to be a constant. It is 
worthwhile to explore some more realistic models in which r varies slowly with time. We 
discuss this possibility in detail here. 

Advocated by the above discussion, one expects that certain amount of interaction can 
alleviate the coincidence problem. One possibility which has been proposed in is to 
allow the ratio of two energy densities to vary slowly but require that there are two positive 
solutions r± to r = 0. Then the coincidence can be alleviated if r_ is close to 0(1) as the 
ratio r evolves from the unstable but finite maximum r. to a stable 
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time, instead of from oo to 0. To demonstrate this possibility we consider again the model 
with an interaction term Q = 36^(1 + r) and a constant u. From equation (fT^. we have 

r' = 3b\l + rf + Siur. (39) 

Obviously setting r' = the equation has two positive solutions with a relation r+r_ = 1. 
It is also possible to show that the ratio will run from an unstable but finite maximum r_|_ 
to a stable minimum r_ at late timejiol. 

However, this does not occur in the context of pure holographic dark energy if one chooses 
the future event horizon as the characteristic size L. In the absence of interaction, the 
dynamics of r is described by 

r' = -r(l + — (40) 

Defining Vl + r = y {y < I) leads to 

2cyV=(l-y')(cy + 2). (41) 

There is only one positive solution to r' = with y = 1 or r = 0. Thus, as in standard 
cosmology r runs from infinity to zero and consequently the coincidence problem still exits 
in this setting. 

Next we intend to propose an alternative way to alleviate the coincidence problem. That 
is, it might not be necessary to have both an unstable finite maximum and a stable minimum 
close to 0(1). The later is more important in the coincidence problem and presumably is 
determined by the physics effective at the current evolution of the universe. The former 
is more related to the early evolution of the universe and whether or not an 0(1) initial 
condition can be obtained is determined by physics beyond the scope of this work. Therefore, 
we would rather leave the question concerning the existence of a positive maximum open 
and concentrate on the models with a positive stable minimum at late time. In particular, 
if we find this stable value is not quite far from the current observation then the coincidence 
problem may be alleviated as the universe has a long time to stay at this stage with a 
similar ratio. Now as an example consider adding the interaction term Q = 36^(1 + r) into 
the holographic dark energy model presented above. Now we find equation (PT|) is changed 
to 

cyy=^{3b'-l)y'-y' + ^y+l. (42) 
11 



Provided c > and 36^ < 1, it still has only one positive solution to y' = but not at r = 0. 
The exact position of the minimum depends on the values of c and b. For explicitness we 
illustrate the evolution of r in FIG. Elwith c = 1 and 6^ = 0.12, which is described by the 
dot-dashed curve. 

Moreover, in the discussion of (|39|) from which two positive solutions are obtained for 
r' = 0, we have made two assumptions, i.e. the coefficients and u are constants. A 
time dependent 6^ or uj might change this situation. In addition, as mentioned in section 
3, specifying any two of L (or pj^), Q, uo and r will determine the other two uniquely. In 
principle, uj might be different from those assumed in the discussion of ()39|1 and consequently 
will lead to different results. Below we will discuss this in more detail through the following 
models. 

A. Model 1: Given Q with time dependent 6^ 

In this model we assume that 

b^ = ble-^/"", 

i.e. 

Q = 362(1 +r)e-'-/^ (43) 

where b"^ =constant. The interaction given by decreases exponentially as r increases 
and consequently at early times when r ^ R, Q is very small and thus can be ignored. 
Therefore, in this model the early age of the universe can be described by the standard 
Friedmann equations without interaction. The interaction becomes important only at late 
time and will lead to a stable minimum which can mitigate the coincidence problem. As 
(fT3|l . one finds that in this case 

r' = 3fe^e-"/^(l + + 3ur, (44) 

and subsequently from (fTHj) one have 

^ = ^[l + u + bl{l + r)e-^/% (45) 

Now to obtain more explicit results we have two options. One is to set to be a constant, 
for example u = —1. This is completely possible in the presence of holographic dark energy 
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and the corresponding size of holographic bound is determined by 

^ = \hl{l + r)e~^-'\ (46) 

With this choice the equation (03} has only one solution to r' = for small R. For example, 
when hi = 0.3, R = 0.25 , the solution to r' = is = 0.196. The second option, instead of 
specifying uj, is setting an appropriate characteristic size L which can also lead to the same 
results obtained above. To show this, we consider the HDE model with L being the future 
horizon. From one has 

3^ = -1 - 3bl{l + r)e-^/« - -=^. (47) 

V 1 + rc 

which lead to 

^ = ^ = -1 + 36^e-^/^i±^ - (48) 

For bl = 0.3, R = 0.5 and c = 1, the solution to r' = is Vf = 0.246. In addition, one can 
check that for r > Vf one always has r < 0. In addition, in this case u also can across — 1 
as in many DE models. For instance, for = 0.1, c = 1 and R = 1, Vf = 0.102, one finds 
that u across — 1 at r = 0.35. 



B. Model 2: Given Q with time independent 

In this subsection we will consider the situation where the dynamics of r at late time can 
be approximated by a power law dependence on a, i.e. 

r = r/ + (ro - r/)a~'', (49) 

under the assumption that Q = 36^(1 + r) where k = 3bl/rf and b^ is a constant. As in 
Section 3, now one can solve for u and p\. From ((7j) and (jH)) it is easy to find that 

u = -blr - (2 + ^)bl. (50) 

Then from ()39|) . one finds that the only solution to r' = is r = rj. Note that this result 

does not dependent on the choice of /c, as what can be obtained from directly. The DE 
density is found to be 

PA = p°a-3(i"^-^^^/). (51) 
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FIG. 1: Distance moduli vs. redshift plot for model 2 (6^ = 0.12, rj = 0.2). The data points are 
from the gold sample of type la supernovae of . 



Then from (|l(jp one finds that 

° 1 + ro 

From (lllj) one finds that the condition a > requires that 

b^> L±I5 . 

' 3(ro + 2 + l/r/) 

As an example, rj = 0.2 leads to 6^ > 0.06. To compare the predictions of this model with 
low redshift observations, the distance moduli vs. redshift are plotted in FIGffl 

As shown in Section 3, given Q and r, the characteristic size L and the dynamics of u 
can be determined uniquely. In fact, from ()51|1 one finds immediately that in HDE models 

L oc J^'-'"-'"/'-f\ (52) 

Moreover, from another point of view, (j49|) and (|50p can also be considered as the conse- 
quences of the characteristic size given by (j52j) . 
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FIG. 2: Distance moduli vs. redshift plot for model 3 (c = 3, A = 0.5). The data points are from 
the gold sample of type la supernovae of [l|. 



C. Model 3: Given the characteristic scale of holographic bound L 



Similar to the previous subsection, here we consider the situation where the late time 
evolution of r to rj can be approximated by an exponential function of a, i.e. r = rj(l + 
7e~'^") where both A and 7 are constants. For simplicity, in the following discussion we set 
7 = 1 and thus we have 

r = r;(l + e-^'^), (53) 

which leads to r/- = ro/(l + e~^\ Nevertheless, rather than fixing the interaction term Q, 
we consider in this subsection the HDE model with the characteristic size L being the future 
event horizon. As discussed in Section 3, from (fT^ and (|^^ . the interaction Q and lo are 
found to be 



Q = r — Xa{r — Vf) + 



2r 



ca/1 + r ' 



1 / 2y/TTr Aae"^" 



(54) 
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From and ()25p one finds that a > and r' < require that 



< l + r + 3cj < 0. 



Then from ()54|) one has 



< — + ^ , _^ ro < 0. 



e 



Since Xae < e ^, 



< . r < e-\ 



1 + e-^ 1 + e-^ 

It is easy to check that the above requirement can be satisfied for any A if c < 2e/ro ~ 15. 
Again we can compare the predictions of this model with low redshift observations, as shown 
in FIGH 

Moreover, to compare the above three models and the interacting HDE model discussed 
at the beginning of this section (see (j42p ). the dynamics of r in these models are plotted in 
FIGElin which the dotted curve corresponds to model 1, the dashed curve to model 2, the 
solid curve to model 3 and the dot-dashed curve to ()42|1 . The parameters used for model 2 
and 3 are the same as those for FIG^and FIGIHI For model 1, the parameters are given in 
the sentence following ()44|). For the curve corresponding to ()42|). c and IP' in ()42|) are taken 
to be 1 and 0.12, respectively. 



VI. DISCUSSION 



We have presented a general formalism of interacting and holographic dark energy in 
this paper. Let us summarize the main results as follows. First we pointed out that in this 
general formalism both the characteristic size of holographic bound L and the coupling term 
of interaction Q for dark energy are not necessarily fixed as in previous references where 
these two sorts of models are separately investigated. Given the conditions that our universe 
is currently accelerating and the ratio of dark matter to dark energy decreases, we derived 
the general relations among the quantities of L,Q,uj and r as well as the constraints on the 
possible range of these quantities. In particular, the dynamics of parameters u and r are 
determined by the choice of L and Q, thus providing a mechanism to change the evolution 
of r from that in standard cosmology such that the coincidence problem may be solved. 
This is the main feature of our formalism. Then we proposed three kinds of strategies to 
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FIG. 3: r vs. log[a] (see text for explanations and parameters used for different curves) 

show how the coincidence problem can be alleviated in this context. One possibility is to 
have a constant ratio throughout the evolution of the universe. The second is to have two 
constant solutions to the ratio r such that it will run from the maximum constant to the 
minimum stable one, while the third and perhaps the most practical one is to have stable 
constant solution at the late time but this value is not quite far from our current observation. 
Focusing on the third strategy we constructed some models explicitly and show how this 
can be implemented by appropriately choosing the quantities of L and Q. In particular our 
results show that at least in the low red-shift region these models are well fitted with the 
current observation. 

In all this paper we assume our universe is spatially flat but it is completely possible 
to show that the parallel analysis could be extended to the spatially closed and hyperbolic 
universe. We also expect that the further investigation will provide us a more exact picture 
of the dark matter and dark energy by strictly fitting the observations in high red-shift 
region. 
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